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TITLE OF THE INVENTION 

FUNCTION ARITHMETIC METHOD AND FUNCTION 
ARITHMETIC CIRCUIT 

5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates generally to 
a function calculation method and a function 
arithmetic circuit for calculating a sine 

10 function and a cosine function with the Taylor 
series equation and, more particularly, to a 
function calculation method and a function 
arithmetic circuit for calculating for use a sine 
function and a cosine function in synchronism with 

15 an arithmetic cycle of fast Fourier transform. 

2. Description of the Related Arts 

Traditionally, in information processing 
apparatuses using computers, sine functions and 
cosine functions calculated from angle 

20 information are sometimes used. Especially, 

when signal processes such as Fourier transform 
are performed by hardware, values of a sine 
function and a cosine function are found using 
angle information as an address and are used in 

25 an arithmetic process. 

Fig. 1 shows an example of a conventional 
system which performs fast Fourier transform in 
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the field of radio astronomy. An analogue 
time-series signal 202 such as electric waves and 
audio signals received by a dedicated parabolic 
antenna 200 are converted by an AD conversion 
5 apparatus into a digital signal 205 and then 
decomposed to a frequency component 208 by 
executing a fast Fourier transform process in a 
Fourier transform apparatus 206 and analyzed. In 
this Fourier transform apparatus 206, angle 

10 information is considered as an address, and a 
sine function and a cosine function corresponding 
to the address are used as arithmetic parameters. 

Figs. 2A and 2B show an example of a 
conventional pipeline fast Fourier transform 

15 apparatus (see Japanese Patent Application 

Lald-Open Pub. No. 06-019955). In Figs. 2A and 
2B, the conventional fast Fourier transform 
apparatus consists of sorting circuit units 242, 
252, 256 and 260, butterfly arithmetic units 244, 

20 254, 258 and 262 and twist coefficient memory 
circuit units 264, 266, 268 and 270, and for 
example, the butterfly arithmetic units 244 
consists of an adder 246 and multipliers 248 and 
250. This fast Fourier transform apparatus is a 

25 circuit performing fast Fourier transform of 

which the number of fast Fourier transform (FFT) 
points to be processed is N, divides an input data 
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point number N into N/R (wherein R Is a radix), 
arranges the divided input data point number in 
serial. Inputs it to a basic circuit of the fast 
Fourier transform which consists of a data sorting 
5 circuit unit of R input, a twist coefficient 

multiplier unit and a butterfly arithmetic unit, 
and is configured such that the fast Fourier 
transform are performed with considering this 
basic circuit as one (1) stage and arranging M 

10 stages (M=logR N) in parallel. In this fast 

Fourier transform, a twist coefficient derived 
from a sine function and a cosine function of angle 
information indicating an address is found in 
advance, and the process is executed with holding 

15 this in memory constituting the twist coefficient 
memory circuit units 264, 266, 268 and 270, 
considering the angle information as an address. 

By the way, traditionally, as a method for 
finding a sine function and a cosine function 

20 corresponding to the angle information 

indicating the address, a look-up table 210 of Fig. 
3 is known. The look-up table 210 consists of 
memory 212 holding values of the sine function and 
the cosine function, considering the angle 

25 information as an address i. In the fast Fourier 
transform apparatus of Figs. 2A and 2B, the twist 
coefficient is calculated from the sine function 



4 



and the cosine function calculated in advance and 
is stored in the memory, but when resolution of 
the angle information becomes higher, capacity of 
memory holding the twist coefficient is increased 
5 in proportion to the resolution. Therefore, if 
the sine function and the cosine function are 
calculated in the fast Fourier arithmetic at a 
high speed, the memory for calculating and storing 
the twist coefficient, and when the resolution is 

10 Improved, capacity of the memory will not 

increased. As these methods for calculating the 
sine function and the cosine function in 
conformity with the arithmetic speed of the fast 
Fourier transform, a circuit in which an 

15 arithmetic process based on a Taylor series 
equation is simply achieved by hardware is 
considered. An arithmetic principle of this 
circuit is as follows. Assuming that angles 0 to 
360 degrees are divided equally into N, the i-th 

20 angle 0 1 is represented by following equation. 
0 1 = 2 7C1/N 

The Taylor series of the sine and the cosine for 
this 1-th angle 0 1 is as follows. 

sin0i = 9i ■ ( 9 i 3 / 3/ ) + ( 0 i 5 I 5/ ) 
25 - ( 0 i 7 / 7/ ) + ( 0 i 9 / 9/ ) - 

= (27C1/N) - {(27T1/N) 3 / 3/ } + {(27T1/N) 5 
/ 5/ } - {(27C1/N) 7 / 7/ } + 
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cos0i = 1 - ( 0 ± 2 / 2/ ) + ( 0 ± 4 / 4/ ) - 
( 0 l 6 / 6/ ) + (0 i 8 / 81 ) - 

= 1 - {(27C1/N) 2 / 2/ > + {(27ri/N) 4 / 41 } 
- {(27C1/N) 6 / 61 > + 
5 Fig. 4 shows a simple circuit based on 

hardware which finds sln0 from the Taylor series, 
and the arithmetic of the Taylor series equation 
up to the fourth term is achieved by a multiplier 
218, a shifter 220, adders 222, 224, 226, 230, 232, 

10 236 and 238 and adders 228, 234, 240. In this 
circuit, for example In the case of the Taylor 
series which has four (4) terms, the value of the 
sine function can be calculated In about 10 clock 
cycles. Further, as a method for quickly 

15 calculating a transcendental function including 
trigonometric functions, there is an apparatus 
which executes parallel processing of the Taylor 
series (see Japanese Patent Application 
Lald-Open Pub. No. 10-214176). This apparatus 

20 executes a parallel process in which the Taylor 
series is divided into two (2) partial series and 
processed in parallel and then two (2) results are 
subtracted. For example. In the case of the 
Taylor series which has 10 terms, by decomposing 

25 it to partial series which have five (5) terms and 
processing them in parallel, the sine function and 
the cosine function can be calculated In about 10 
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clock cycles. Each parallel processing unit 
executes the process by transforming cyclically, 
and floating-point arithmetic is assumed. 

However, in these conventional methods for 
5 finding a sine function and a cosine function with 
the use of Taylor series equations, there are 
following problems. First, as shown in Fig. 4, 
in the case that each term of the Taylor series 
equation is simply achieved by a hardware 

10 arithmetic circuit, scale of hardware is 

Increased corresponding to the number of terms, 
and clock cycles required for the arithmetic are 
also increased, therefore there is a problem which 
reduce the overall arithmetic speed if it is 

15 applied to fast Fourier transform. On the other 
hand, in the method which calculates in parallel 
by dividing the Taylor series into two (2) partial 
series, although the arithmetic process speed can 
be Increased by executing it in parallel, the 

20 scale of hardware is equivalent to the case that 
it is not divided into partial series, and the 
scale of circuit is large, and if achieved by 
software, there is a problem of a high processing 
load, because parallel processing of 

2 5 floating-point arithmetic is executed. 
SUMMARY OF THE INVENTION 

According to the present invention there are 
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provided a function arithmetic method and a 
function arithmetic circuit for calculating a 
sine function and a cosine function at a high speed 
utilizing the Taylor series equation. The 
5 present invention aims at providing a function 
arithmetic method and apparatus for calculating 
a sine function and a cosine function at a high 
speed utilizing the Taylor series equation, 
without using look-up table memory, for use in 
10 calculation of a twist coefficient of fast Fourier 
transform . 

The present invention provides a function 
arithmetic method for calculating a sine function 
sin 0 1 using a Taylor series equation. The 

15 function arithmetic method comprises a cyclic 
equation setting step, with a cyclic equation 
setting unit, transforming and setting a Taylor 
series equation for calculating a sine function 
Into a single cyclic equation (Q ■ K+S'X 2 *Q) common 

20 to terms of the Taylor series equation, the single 
cyclic equation having a new known number Q that 
Is defined by multiplying a known number Q and the 
square of a variable X, shifting the result by a 
shift number S and then adding a constant K 

25 thereto; an adjustment step, with an adjustment 
unit, adjusting and preparing the shift number S 
such that within a variation range of the variable 



X the variable X has a maximum value 1 with the 
constant K being not greater than 1; and a cyclic 
equation executing step, with a cyclic equation 
execution unit, inputting and converting angle 
5 information 1 to the variable X, and executing the 
cyclic equation in sequence from higher order term 
to lower order term for the number of terms of the 
Taylor series equation to derive a sine function 
of the angle information 1. In this case, the 

10 cyclic equation executing step includes 

executing the arithmetic process of the known 
number Q, the variable X, an intermediate value 
after the shifting and the constant K of the cyclic 
equation with the number of bits that is obtained 

15 by adding the number of protect bits to the final 
number of the bits* 

The present invention provides a function 
arithmetic method for calculating a cosine 
function cos d 1 using a Taylor series equation. 

20 The function arithmetic method comprises a cyclic 
equation setting step, with a cyclic equation 
setting unit, transforming and setting a Taylor 
series equation for calculating a cosine function 
into a single cyclic equation (Q - K+S-X 2 'Q) common 

25 to terms of the Taylor series equation, the single 
cyclic equation having a new known number Q that 
is defined by multiplying a known number Q and the 
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square of a variable X, shifting the result by a 
shift number S and then adding a constant K 
thereto; an adjustment step, with an adjustment 
unit, adjusting and preparing the shift number 
5 S such that within a variation range of the 

variable X the variable X has a maximum value 1 
with the constant K being not greater than 1; and 
a cyclic equation executing step, with a cyclic 
equation execution unit, inputting and 

10 converting angle information i to the variable X, 
and executing the cyclic equation in sequence from 
higher order term to lower order term for the 
number of terms of the Taylor series equation to 
derive a cosine function of the angle information 

15 i. In this case as well, the cyclic equation 

executing step includes executing the arithmetic 
process of the known number Q, the variable X, an 
intermediate value after the shifting and the 
constant K of the cyclic equation with the number 

20 of bits that is obtained by adding the number of 
protect bits to the final number of the bits. 

The present invention provides a function 
arithmetic circuit for calculating a sine 
function sine 0 1 using a Taylor series equation. 

25 The function arithmetic circuit comprises a 

cyclic equation arithmetic unit calculating a 
cyclic equation (Q = K+S • X 2 • Q) that is obtained 
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by transforming a Taylor series equation for 
calculating a sine function, the cyclic equation 
having a new known number Q that is defined by 
multiplying a known number Q and the square of a 
variable X, shifting the result by a shift number 
S and then adding a constant K thereto; a 
conversion adjustment unit converting input 
angle Information 1 into the variable X, as well 
as adjusting and outputtlng the shift number S 
such that the variable X has a maximum value 1 
within a variation range of the variable X; a 
constant table finding in advance and holding 
constants K corresponding to terms of a Taylor 
series equation for calculating a sine function 
and the shift numbers adjusted such that the 
constants K becomes not greater than 1; and an 
arithmetic control unit causing the cyclic 
equation arithmetic unit to perform a cyclic 
arithmetic in sequence, based on the selection of 
the constant K and the shift number S of the 
constant table, from higher order term to lower 
order term for the number of terms of the Taylor 
series equation defined In advance when the 
variable X is output from the conversion 
adjustment unit, to thereby derive a sine function 
of the angle Information 1. In this case, the 
cyclic equation arithmetic unit executes the 
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arithmetic process of the known number Q, the 
variable X, an intermediate value after the 
shifting and the constant K of the cyclic equation 
with the number of bits that is obtained by adding 
5 the number of protect bits to the final number of 
the bits. 

The present invention provides a function 
arithmetic circuit for calculating a cosine 
function cos 0 i using a Taylor series equation. 

10 The function arithmetic circuit comprises a 

cyclic equation arithmetic unit calculating a 
cyclic equation (Q = K+S " X 2 • Q) that is obtained 
by transforming a Taylor series equation for 
calculating a cosine function, the cyclic 

15 equation having a new known number Q that is 

defined by multiplying a known number Q and the 
square of a variable X, shifting the result by a 
shift number S and then adding a constant K 
thereto; a conversion adjustment unit converting 

20 input angle information 1 Into the variable X, as 
well as adjusting and outputting the shift number 
S such that the variable X has a maximum value 1 
within a variation range of the variable X; a 
constant table finding in advance and holding 

25 constants K corresponding to terms of the Taylor 
series equation for calculating a cosine function 
and the shift numbers adjusted such that the 
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constants K become not greater than 1; and an 
arithmetic control unit causing the cyclic 
equation arithmetic unit to perform a cyclic 
arithmetic in sequence, based on the selection of 
5 the constant K and the shift number S of the 

constant table, from higher order term to lower 
order term for the number of terms of the Taylor 
series equation defined in advance when the 
variable X is output from the conversion 

10 adjustment unit, to thereby derive a cosine 

function of the angle information 1. In this case, 
the cyclic equation arithmetic unit executes the 
arithmetic process of the known number Q, the 
variable X, an intermediate value after the 

15 shifting and the constant K of the cyclic equation 
with the number of bits that is obtained by adding 
the number of protect bits to the final number of 
the bits. 

The present invention provides a function 
20 arithmetic circuit for pipeline calculating a 
sine function sin 0 i using a Taylor series 
equation. The function arithmetic circuit 
comprises a pipeline arithmetic unit forming a 
pipeline connection which Includes cyclic 
25 equation arithmetic units each provided for each 
term and calculating a cyclic equation (Q = K+S • 
X 2 • Q) obtained by transforming a Taylor series 
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equation for calculating a sine function, the 
cyclic equation having a new known number Q that 
is defined by multiplying a known number Q and the 
square of a variable X, shifting the result by a 
5 shift number S and then adding a constant K 

thereto; a conversion adjustment unit converting 
input angle information i into the variable X and 
adjusting the shift number S such that the 
variable X has a maximum value 1 within a variation 

10 range of the variable X for the output to the 
pipeline arithmetic unit; a constant table 
finding in advance and holding the constants K 
corresponding to terms of the Taylor series 
equation for calculating a sine function and the 

15 shift numbers adjusted such that the constants K 
become not greater than 1; and a pipeline control 
unit causing the cyclic equation arithmetic units 
of the pipeline arithmetic unit to select the 
constant K and the shift number S of the 

20 corresponding term of the Taylor series equation 
from the constant table, to calculate in parallel 
and to derive a sine function of the angle 
information i based on the output of the cyclic 
equation arithmetic unit at the final stage, each 

25 time the variable X is output from the conversion 
adjustment unit. In this case, the cyclic 
equation arithmetic units of the pipeline 
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arithmetic unit execute the arithmetic process of 
the known number Q, the variable X, an 
intermediate value after the shifting and the 
constant K of the cyclic equation with the number 
of bits that Is obtained by adding the number of 
protect bits to the final number of the bits. 

The present invention provides a function 
arithmetic circuit for pipeline calculating a 
cosine function cos 6 i using a Taylor series 
equation. The function arithmetic circuit 
comprises a pipeline arithmetic unit forming a 
pipeline connection which includes cyclic 
equation arithmetic units each provided for each 
term and calculating a cyclic equation (Q = K+S • 
X 2 " Q) obtained by transforming a Taylor series 
equation for calculating a cosine function, the 
cyclic equation having a new known number Q that 
is defined by multiplying a known number Q and the 
square of a variable X, shifting the result by a 
shift number S and then adding a constant K 
thereto; a conversion adjustment unit converting 
input angle information 1 into the variable X and 
adjusting the shift number S such that the 
variable X has a maximum value 1 within a variation 
range of the variable X for the output to the 
pipeline arithmetic unit; a constant table 
finding in advance and holding the constants K 
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corresponding to terms of the Taylor series 
equation for calculating a cosine function and the 
shift numbers adjusted such that the constants K 
become not greater than 1; and a pipeline control 
5 unit causing the cyclic equation arithmetic units 
of the pipeline arithmetic unit to select the 
constant K and the shift number S of the 
corresponding term of the Taylor series equation 
from the constant table, to calculate in parallel 

10 and to derive a cosine function of the angle 

information i based on the output of the cyclic 
equation arithmetic unit at the final stage, each 
time the variable X is output from the conversion 
adjustment unit. In this case as well, the cyclic 

15 equation arithmetic units of the pipeline 

arithmetic unit execute the arithmetic process of 
the known number Q, the variable X, an 
intermediate value after the shifting and the 
constant K of the cyclic equation with the number 

20 of bits that Is obtained by adding the number of 
protect bits to the final number of the bits. 

The present invention provides a function 
arithmetic circuit for pipeline calculating a 
sine function and a cosine function, suitable for 

25 a fast Fourier transform apparatus. The function 
arithmetic circuit comprises a sine function 
arithmetic circuit and a cosine function 
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arithmetic circuit. The sine function 
arithmetic circuit includes a pipeline 
arithmetic unit forming a pipeline connection 
which Includes cyclic equation arithmetic units 
5 each provided for each term and calculating a 
cyclic equation (Q = K+S • X 2 • Q) obtained by 
transforming a Taylor series equation for 
calculating a sine function, the cyclic equation 
having a new known number Q that is defined by 

10 multiplying a known number Q and the square of a 
variable X, shifting the result by a shift number 
S and then adding a constant K thereto; a 
conversion adjustment unit converting input 
angle information 1 into the variable X and 

15 adjusting the shift number S such that the 

variable X has a maximum value 1 within a variation 
range of the variable X for the output to the 
pipeline arithmetic unit; a constant table 
finding in advance and holding the constants K 

20 corresponding to terms of the Taylor series 

equation for calculating a sine function and the 
shift numbers adjusted such that the constants K 
become not greater than 1; and a pipeline control 
unit causing the cyclic equation arithmetic units 

25 of the pipeline arithmetic unit to select the 
constant K and the shift number S of the 
corresponding term of the Taylor series equation 
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from the constant table, to calculate ±n parallel 
and to derive a sine function of the angle 
Information ± based on the output of the cyclic 
equation arithmetic unit at the final stage, each 
time the variable X is output from the conversion 
adjustment unit. The cosine function arithmetic 
circuit includes a pipeline arithmetic unit 
forming a pipeline connection which Includes 
cyclic equation arithmetic units each provided 
for each term and calculating a cyclic equation 
(Q = K+S • X 2 ■ Q) obtained by transforming a Taylor 
series equation for calculating a cosine function, 
the cyclic equation having a new known number Q 
that is defined by multiplying a known number Q 
and the square of a variable X, shifting the result 
by a shift number S and then adding a constant K 
thereto; a conversion adjustment unit 
converting input angle information i into the 
variable X and adjusting the shift number S such 
that the variable X has a maximum value 1 within 
a variation range of the variable X for the output 
to the pipeline arithmetic unit; a constant table 
finding in advance and holding the constants K 
corresponding to terms of the Taylor series 
equation for calculating a cosine function and the 
shift numbers adjusted such that the constants K 
become not greater than 1; and a pipeline control 
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unit causing the cyclic equation arithmetic units 
of the pipeline arithmetic unit to select the 
constant K and the shift number S of the 
corresponding term of the Taylor series equation 
5 from the constant table, to calculate in parallel 
and to derive a cosine function of the angle 
Information 1 based on the output of the cyclic 
equation arithmetic unit at the final stage, each 
time the variable X is output from the conversion 

10 adjustment unit. In this case, the cyclic 
equation arithmetic units of the pipeline 
arithmetic unit execute the arithmetic process of 
the known number Q, the variable X, an 
intermediate value after the shifting and the 

15 constant K of the cyclic equation with the number 
of bits that Is obtained by adding the number of 
protect bits to the final number of the bits. 

In this function arithmetic circuit, twist 
coefficient values of a plurality of butterfly 

20 stages provided in a pipeline fast Fourier 
transform apparatus whose radix Is 2 are 
calculated based on the sine function and cosine 
function of the input information 1. 

The present invention provides a function 

25 arithmetic method for calculating, using a Taylor 
series equation , general transcendental 
functions including sine functions and cosine 
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functions, I.e., transcendental functions that 
are defined as functions that cannot be expressed 
In the form of polynomials, such as trigonometric 
functions, logarithms and exponential functions. 
5 The function arithmetic method comprises a cyclic 
equation setting step, with a cyclic equation 
setting unit, transforming and setting a Taylor 
series equation for calculating a transcendental 
function into a single cyclic equation (Q = K+S ■ 

10 X 2 *Q) common to terms of the Taylor series equation, 
the cyclic equation having a new known number Q 
that is defined by multiplying a known number Q 
and a variable X, shifting the result by a shift 
number S and then adding a constant K thereto; an 

15 adjustment step, with an adjustment unit, 

adjusting and preparing the shift number S such 
that within a variation range of the variable X 
the variable X has a maximum value 1 with the 
constant K being not greater than 1; and a cyclic 

20 equation executing step, with a cyclic equation 
execution unit, converting Input information to 
the variable X and executing the cyclic equation 
in sequence from higher order term to lower order 
term for the number of terms of the Taylor series 

25 equation to thereby derive a transcendental 
function of the input information. 

Thus, according to the present invention, by 
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transforming a Taylor series equation into one (1) 
cyclic equation which is common to terms of the 
Taylor series equation, details of arithmetic 
which is performed from higher order to lower 
5 order can be standardized, and in both cases of 
performing cyclic calculation and performing 
pipeline calculation, an arithmetic algorithm is 
simplified; a processing load is reduced in 
software; a circuit structure is simplified in 

10 hardware; and a higher speed can be achieved in 
both. Also, by introducing shift operations of 
bits, even with the same minimum small number of 
arithmetic bits, necessary arithmetic bit 
accuracy can be ensured. As a result, without 

15 using floating-point arithmetic or an arithmetic 
circuit with the large number of bits, fixed-point 
arithmetic circuit with the small number of bits 
can be used, and function calculation with small 
hardware can be easily achieved. Further, by 

20 applying a pipeline function arithmetic circuit 
of the present invention to a fast Fourier 
transform apparatus, a twist coefficient can be 
calculated in real time, and conventional memory 
for holding the twist coefficient is unnecessary, 

25 and if the number of fast Fourier transform points 
becomes larger, the twist coefficient can be 
accurately generated by smaller hardware than the 
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case of using the twist coefficient memory. 

Consequently, although the conventional 
Fourier transform apparatus is constructed by 
Integrating butterfly arithmetic circuits and 
5 others into LSI and using external memory for the 
twist coefficient, it is possible to relatively 
easily construct it with LSI only, by using the 
pipeline function arithmetic circuit of the 
present invention . 

10 The above and other objects, features and 

advantages of the present invention will become 
more apparent from the following detailed 
description with reference to the drawings. 
BRIEF DESCRIPTION OF THE DRAWINGS 

15 Fig. 1 is an explanatory view of a system 

using a conventional fast Fourier transform 
apparatus ; 

Figs. 2A and 2B are block diagrams of a 
conventional fast Fourier transform apparatus 
20 holding a twist coefficient in memory; 

Fig. 3 is an explanatory view of a look-up 
table which converts angle Information address i 
Into a sine function or a cosine function; 

Fig. 4 is a circuit diagram as simple hardware 
25 which calculates a sine function by the use of the 
Taylor series equation; 

Fig. 5 is a block diagram of a functional 
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structure of a function arithmetic process of the 
present invention achieved by a program of a 
computer; 

Fig. 6 is an explanatory view of table 
5 contents used in function arithmetic of Fig. 5; 

Fig. 7 is a flowchart of a function arithmetic 
process of Fig. 5; 

Fig. 8 is a flowchart of a sine function 
arithmetic process of Fig. 5; 
10 Fig. 9 is a flowchart of a cosine function 

arithmetic process of Fig. 5; 

Fig. 10 Is a block diagram of a function 
arithmetic circuit of the present invention 
achieved by hardware; 
15 Fig. 11 is a circuit block diagram of a sine 

function arithmetic circuit provided in Fig. 10; 

Fig. 12 is a circuit block diagram of a cosine 
function arithmetic circuit provided in Fig. 10; 

Fig. 13 Is a flowchart of a sine function 
20 arithmetic process according to the present 
invention which keeps accuracy by rounding 
arithmetic results to which protect bits are 
added; 

Fig. 14 is a flowchart of a cosine function 
25 arithmetic process according to the present 
invention which keeps accuracy by rounding 
arithmetic results to which protect bits are 
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added; 

Fig. 15 Is a circuit block diagram of a sine 
function arithmetic circuit of the present 
invention which performs pipeline arithmetic ; 
5 Fig. 16 is a circuit block diagram of a cosine 

function arithmetic circuit of the present 
invention which performs pipeline arithmetic; 

Figs. 17A and 17B are block diagrams of a fast 
Fourier transform circuit whose radix is 2 , to 
10 which a function arithmetic circuit of the present 
invention performing the pipeline arithmetic of 
Fig. 15 and Fig. 16 as twist coefficient 
arithmetic is applied; 

Fig. 18 is an explanatory view of sorting 
15 across four (4) stages of Figs. 17A and 17B and 
I/O data of butterfly arithmetic ; 

Figs. 19A and 19B are block diagrams in which 
the fast Fourier transform circuit of Figs. 17A 
and 17B is in a parallel configuration; and 
20 Fig. 20 Is a block diagram of a twist 

coefficient arithmetic circuit used in the fast 
Fourier transform circuit of Figs. 19A and 19B. 
DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 5 is a block diagram of a functional 
25 structure of a function arithmetic process of the 
present Invention achieved by a program of a 
computer. In Fig. 5, a function arithmetic 
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processing system of the present Invention 
consists of an angle detection unit 10 and a 
computer 11 As the angle detection unit 10, for 
example, a rotary encoder and others are used, 
5 which output address i when angles 0 to 360 degrees 
are divided equally into N and represent an angle 
6 i of this address 1. In the computer 11, a 
function arithmetic processing unit 12 is 
provided, which Inputs the address i from the 

10 angle detection unit 10 and calculates a sine 
function and cosine function of sin 0 1 and cos 
0 1. The function arithmetic processing unit 12 
consists of a cyclic equation setting unit 13, an 
adjustment unit 14 and a cyclic equation executing 

15 unit 15, and a constant table 16 is further 

provided, which stores constant and others used 
by the cyclic equation executing unit 15. The 
cyclic equation setting unit 13 transforms a 
Taylor series equation for calculating the sine 

20 function and the cosine function into a single 
cyclic equation which is common to terms of the 
Taylor series equation of the Taylor series 
equation and set this to the cyclic equation 
executing unit 15. This cyclic equation is given 

25 by following equation. 

Qd = Kd + S(d+2) -Xi 2 - [Q(d+2) ] (1) 
In the cyclic equation (1), the Taylor series 
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equation is to multiply known number Q(d+2) by a 
squared variable Xi and shift the result by a shift 
number S(d+2). A transform process from this 
Taylor series equation to the cyclic equation is 
5 clarified in description below. For the cyclic 
equation of equation (1), the adjustment unit 14 
adjusts the variable XI such that a maximum value 
of the variable XI becomes 1 within a transform 
range and the shift number S(d+2) such that the 

10 a constant Kd will be not greater than 1 and 

prepares this in the constant table 16 in advance. 
In the case that the sine function and the cosine 
function are calculated for angles of 0 to 45 
degrees and the case that calculated for angles 

15 of 0 to 90 degrees. An example of the constant Kd 
and the shift number Sd stored in the constant 
table 16 are shown in Fig. 6. 

In the constant table 16 of Fig. 6, values 
of the constant Kd and the shift numbers Sd in the 

20 cyclic equation of equation (1) are shown 

corresponding to the parameters d and (27t) d /d! 
of the Taylor series equation for each of the sine 
function and the cosine function, and the shift 
numbers Sd are divided into the case of 

25 calculating up to 45 degrees and the case of 
calculating up to 90 degrees and stored. 

Fig. 7 is a flowchart of a process procedure 
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of the function arithmetic processing unit 12 
according to the present invention achieved by- 
program control of the computer 11 of Fig. 5. 

In this function arithmetic process, after 
5 a setting process of the cyclic equation of the 
equation (1) is performed in step SI, by adjusting 
the shift number Sd in order to optimize the 
variable XI and the constant Kd in step 2 , the 
table 16 such as shown In Fig. 6 are generated. 

10 If preparation processes in steps SI and S2 are 
completed, an Input of the angle Information 
address 1 from the angle detection unit 10 Is 
checked in step S3, and If there is the Input of 
the address 1, proceed to step S4, and the address 

15 1 as the angle Information is converted to the 
variable XI, and the sine function slnd 1 for the 
angle 6 1 of the address 1 is calculated by 
executing the cyclic equation from higher order 
term to lower order term for the number of terms 

20 of the Taylor series equation. A transform 

process from the Taylor series equation to the 
cyclic equation of the equation (1) is then 
described. In the function arithmetic process of 
the present invention is divided into the case 

25 that the sine function and the cosine function are 
calculated up to an angle of 4 5 degrees and the 
case that these are calculated up to an angle of 
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90 degrees. At this point, the sine function and 
the cosine function in the range of angles 90 to 
180 degrees can be easily found from the sine 
function of angles 0 to 90 degrees which is the 
5 arithmetic result, and the function arithmetic 
according to the present invention is not needed 
for these. In other words, formulas of 
trigonometric functions, which are: 

cos(a + /3) = cosa -cos/3 - sina -sin/3 
10 sin(a + /3) a cosa a sinj3 + sina -cos/3 

can be used to derive these. These formulas of 
trigonometric functions are transformed into 
following equations. 

oos(90° + /3)=cos90° *cos/3 - sin90° -sin/S = 
15 0'oos/3 - l*sln/3 = - sin/3 

sln(90° + /3)=oos90° *sin/3 + sin90° 'cos/3 = 
O'sinjS + l-cos/3 =a cos/3 

For example, in the case of an angle of 120 degrees, 
since 90° +/3=120° , /3 =30° is derived. 

20 Therefore: 

cosl20° =-sin30° ; and 
slnl20° =oos30° 
can be found. The sine function and the cosine 
function in the range of angles 45 to 90 degrees 

25 can also be easily found from the sine function 
and the cosine function of angles 0 to 45 degrees 
which is found by the function arithmetic of the 
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present Invention. In other words, these can be 
derived from following equations which are 
formulas of trigonometric functions. 

cos(CK-j3) = cos<Z*cosj3 + slntt*slnj3 
5 sln(a-j3) = -cosa-sinjS + sinCK •cosj3 

These formulas of trigonometric functions are 
transformed into following equations. 

cos(90° - /3 ) = cos90° # cosj3 + sin90° -sin/3 = 
0*cos/3 + l*slnj3 = sinjS 
10 sin(90° - j3 ) = -cos90° -sinj9 + sin90° -cos 

jB « -0*sinj3 + l*cosj3 = cos/3 

For example. In the case that the sine function 
and the cosine function of an angle of 60 degree, 
since 90° - j3 =60° , 3 =30° is derived. 

15 Therefore z 

cos60° =sin30 ° ; and 
sin60° =cos30° 
can be found. 

Described here is transformation to the 

20 cyclic equation in the case that the sine function 
and the cosine function for angles 0 to 45 degrees 
are found. Since the angle address 1 is within 
the range of angles 0 to 45 degrees, the value of 
the address varies from 0 to N/8. Since the 

25 address 1 represents the 1-th angle 0 1 in the case 
that angles 0 to 360 degrees are divided equally 
into N, this is given by following equation. 
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6 ± = 27T1/N (2) 
Using the equation (2) which indicates the angle 

0 1, the Taylor series equation is developed as 
follows . 

5 sln0i =01 - ( 0 i 3 / 3/ ) + ( 0 i 5 / 51 ) 

- ( 0 i 7 / 7/ ) + ( 0 i 9 / Si ) - 

= (27C1/N) - {(27C1/N) 3 / 3/ > + {(27C1/N) 5 
/ 5/ } - {(27C1/N) 7 / 11 } + 

- 27ti/{8' (N/8) } - [27Ti/{8- (N/8) }] 3 / 3/ + 
10 [27Ci/{8- (N/8) }] 5 / 5/ 

- {27Ci/{8-(N/8)}] 7 / 7/ + 

- C±/(N/8)}-27T/ 2 3 - {l/(N/8)} 3 -{(27C ) 3 / 
3/ >/ 2 9 + {i/(N/8)} 5 -{(27T ) 5 / 5/ }/ 2 15 

- {i/(N/8)} 7 -{(27T) 7 / 11 }/ 2 21 + 

15 In this development of the Taylor series equation, 
a first row is the original Taylor series equation, 
and in a second row, the value of 0 1 of the 
equation (2) is substituted in the equation. In 
third and fourth rows of the equation, by 

20 substituting N=8(N/8), it is transformed into 
exponentiation of 2 . In fifth and sixth rows, 
portions in brackets are divided into the variable 
portions and the constant portions. 

In the present invention, the constant Kd of 

25 the cyclic equation of the equation (1) will be 
a value which is not greater than 1 and close to 

1 as much as possible. For example, by dividing 
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2 7C=6.28 ±n the first term of fifth and sixth rows 
in the transformed equation of the Taylor series 
equation by 2 3 = 8 , K1 = 0 .785 can be derived. In 
other words, if 2 7C=6.28 is divided by 2 2 = 4 # 
Kl=1.57 is derived, which can not be adopted 
because this is greater than 1, and if divided by 
2 4 =16, Kl=0.3925 is derived, which is far from 1, 
and therefore, by dividing by 2 3 =8, the constant 
Kl can be a value which is not greater than 1 and 
close to 1 as much as possible. The constants of 
the transformed equation of the Taylor series 
equation can be summarized as following equation. 

sinfli = {i/(N/8) }• [27C/ 2 3 ] - { 1/ (N/8 ) } 3 • [ ( 2 
7C) 3 /{3! -2 6 )]/ 2 3 

+ {i/(N/8)} 5 -[(27T ) 5 /{Si -2 7 )]/ 2 8 - 
{i/(N/8)} ? - [ (27C ) 7 /(7/ -2 7 )]/ 2 14 + 
For this equation, the constants Rl, K3 , K5 and 
K7 are defined as follows. 

Kl = [27C/ 2 3 ] 

K3 - - [(27T) 3 /(3/ -2 6 )] 

K5 = [ (27T ) 5 /(5/ -2 ? ) ] 

K7 - - [(27T) 7 /{It -2 7 )] 
By transforming, using the constants Kl, K3 , K5 
and K7 , following equation can be obtained. 

sin0i - {i/(N/8)}-Kl + { 1/ (N/8 ) } 3 -K3/ 2 3 + 
{i/(N/8) } 5 -K5/ 2 8 + {i/(N/8) } ? -K7/ 2 14 + 

- {l/(N/8)}-[Kl + Ci/(N/8) } 2 *K3/ 2 3 + 
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{±/(N/8)} 4 *K5/ 2 8 + {l/(N/8) } 6 'K7/ 2 14 + — ] 
= {i/(N/8)}-[Kl + 2" 3 -{i/(N/8) } 2 ' [K3 + 

{1/(N/8)> 2 'K5/ 2 5 

+ {i/(N/8) } 4 'K7/ 2 11 + — ] ] 
5 = {±/(N/8) }• [Kl + 2 _3 '{±/(N/8) } 2 * [K3 + 2~ 5 ' 

{±/(N/8)> 2 - [K5 

+ {l/(N/8) > 2 'K7/ 2 6 + — ] ] ] 

o £i/(N/8) }* [Kl + 2 _3 -{l/(N/8) > 2 ' [K3 + 2~ 5 • 
{l/(N/8)> 2 ' [K5 

10 + 2 _6 '{l/(N/8) > 2 * [K7 +—]]]] 

In this transformed equation, the constants Kl, 
K3, K5 and K7 are substituted In a first row, and 
{l/(N/8)> are bundled out for the portion in 
brackets In a second row. In third and fourth 

15 rows, 2' 3 Is bundled for first and second terms 
In braces on the right side. In fifth and sixth 
rows, 2' 5 Is further bundled for the second term 
and subsequent terms In that. Finally, In a 
seventh row, 2" 6 is bundled for the third term. 

20 For the above transformed equation, when as 
variable XI: 

XI = l/(N/8) 
Is defined, and as bit shift numbers S3, S5 and 
S7 : 



S3 = 


2~ 


3 


S5 = 


2" 


5 


S7 = 


2" 


6 
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are defined, and as known numbers Q7 , Q5 , Q3 and 
Ql: 

Q7 = K7 + 

Q5 = K5 + S7-X1 2 -[Q7] 
5 Q3 = K3 + S5-Xi 2 *[Q5] 

Ql » Rl + S3-X± 2 -[Q3] 
are defined, following equation can be obtained. 

sin01 = Xi-[R1 + S3-Xi 2 -[K3 + S5-Xi 2 -[K5 + 
S7-X1 2 - [R7 +—]]]] 
10 = Xi-[R1 + S3-Xi 2 -[R3 + S5'X1 2 -[R5 + S7-X1 2 - 

[Q7] ] ] ] 

- Xi*[Rl + S3'X1 2 -[R3 + S5'X1 2 -[Q5] ] ] 
= Xi-[R1 + S3-Xi 2 -[Q3] ] 

- XI- [Ql] 

15 The cyclic portion of the sine function In this 
transformed equation takes values of d= 1, 3, 5, 
7 and following equation is given, which is 

the same as the equation (1). 

Qd » Kd + S(d+2) -XI 2 - [Q(d+2) ] 

20 Therefore, In calculation of slnd unit with the 
use of this cyclic equation, the sine function sin 
0 i can be found by sequentially calculating the 
cyclic equation in the order of d=7, 5, 3 and 1 
to find Q7, Q5 , Q3 and Ql and finally multiplying 

25 by Xi. For Q7 as an initial value, the constant 
K7 Is used. 

Similarly, transform from the Taylor series 
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equation to the cyclic equation for the cosine 
function is described as follows. By- 
substituting the equation (2) in the Taylor series 
equation of cos 0 1, the transformed equation 
becomes as follows. 

cosfll = 1 - ( 6 l 2 / 2/ ) + ( 6 i 4 / 41 ) - 
( 6 i 6 / 61 ) + ( 9 l 8 / 81 ) - 

- 1 - {(27C1/N) 2 / 2/ } + {(27C1/N) 4 / 41 } 
- {(27C1/N) 6 / 61 } + 

- 1 - [27Ci/{8- (N/8)}] 2 / 21 + [27Ci/{8- 
(N/8)}] 4 / 41 

- {27Ci/{8*(N/8)}] 6 / 61 + 

= 1 - {l/(N/8) } 2 -{(27C ) 2 / 21 }/ 2 6 + 
{i/(N/8)} 4 '£(27C ) 4 / 41 }/ 2 12 

- {i/(N/8)} 6 '{(27C ) 6 / 61 }/ 2 18 + 

= 1 - {i/(N/8)} 2 '[(27C) 2 /(21 *2 5 )]/ 2 1 + 
(i/(N/8)} 4 ' [ (27C ) 4 /(41 -2 7 )]/ 2 5 

- {l/(N/8)} 6 « [(2^ ) 6 /(61 -2 7 )]/ 2 11 + 

By defining constants KO, K2 , K4 and K6 as: 

KO - 1 

K2 = - [(27T) 2 /(21 -2 5 )] 

K4 = [(27C) 4 /(41 *2 7 )] 

K6 = - [ (2^ ) 6 /(61 -2 7 ) ] 
and transforming, following equation is given. 

cos 0 i = KO + {i/(N/8) } 2 *K2/ 2 1 + (i/(N/8) } 4 ' 
K4/ 2 5 + {l/(N/8) }*-K6/ 2 11 + 

- KO + 2~ 1 *{±/(N/8) } 2 * [K2 + {i/(N/8) } 2 'K4/ 2 4 
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+ {i/(N/8) } 4 *K6/ 2 10 + — ] 

- KO + 2" 1 ' {i/(N/8) > 2 ' [K2 + 2~ * • { 1/ ( N/8 ) } 2 • 
[K4 + {±/(N/8) } 2 *K6/ 2 6 + — ] ] 

= KO + 2 _1 '{±/(N/8) } 2 ' [K2 + 2 _4 -{l/(N/8)} 2 ' 
5 [K4 + 2" 6 -{i/(N/8) } 2 * [K6 + — ] ] ] 

At this point, by defining the variable Xi as: 

Xi = i/(N/8) 
and the shift numbers S2, S4 and S6 as: 
S2 = 2 _1 
10 S4 = 2~ 4 

S6 = 2" 6 

and the known number Q6 , Q4 , Q2 and Q0 as: 
Q6 = K6 + — 
Q4 ■ K4 + S6'Xi 2 *[Q6] 
15 Q2 ■ K2 + S4'X1 2 '[Q4] 

Q0 = KO + S2'X1 2 '[Q2] 
it is possible to transform to following equation. 

oos0 1 = KO + S2*X1 2 '[K2 + S4'Xi 2 *[K4 + S6" 
XI 2 ' [K6 + — ] ] ] 
20 - KO + S2-X1 2 ' [K2 + S4*Xi 2 ' [K4 + S6'Xi 2 * [Q6] ] ] 

- KO + S2'Xi 2 *[K2 + S4*Xi 2 -[Q4] ] 

- KO + S2'Xi 2 ' [Q2] 
= Q0 

The cyclic portion of the cosine function in this 
25 transformed equation of the Taylor series 

equation takes values of d= 0, 2, 4, 6 *'**, and 
the cyclic equation which is common to terms of 
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the Taylor series equation is: 

Qd = Kd + S(d+2) "Xi 2 - [Q(d+2) ] 
which is the same as the equation (1). 

Fig. 8 is a flowchart of a sine function 
5 arithmetic process in the cyclic equation 

executing unit 15 of Fig, 5, and the process 
procedure is as follows. 

Step Sis Input the address 1 which indicates 
angles . 

10 Step S2: Since calculation from 0 to 

45degrees is performed with the address 1, define 
the variable Xi as Xi=i/(N/8) and normalize such 
that the maximum value results in 1. 

Step S3: Calculate the square of the variable 

15 Xi. 

Step S4: Set d=7 which is the number of terms 
Indicating termination of the Taylor series 
equation, find the constant K(7) from the constant 
table 16 of Fig. 6 and substitute this for the 
2 0 known number Q ( 7 ) . 

Step S5 : Define d as d=d-2. 

Step S6: Check whether d is not less than 1. 
If it is not less than 1, proceed to step S7„ and 
if less than 1, proceed to step S8. 
25 Step S7: Sequentially calculate the cyclic 

equation in the order of D=5, 3, 1, through the 
route of steps S5 to S7, for the known number Q5 , 



Q3 and Ql . 

Step S8: Tills Is the case that arithmetic of 
the cyclic equation Is terminated because d 
becomes less than 1, and therefore , by multiplying 
5 the known number Q(l) by the variable XI, sin 

0 1 Is calculated. 

Step S9 : Output the calculated value of sin 

01. 

Fig. 9 is a flowchart of a cosine function 
10 arithmetic process of the present invention of Fig. 
5. The function arithmetic process is as 
follows. 

Step SI: Input the address 1 which indicates 
angles. 

15 Step S2: Divide the address i by (N/8) to find 

the variable XI which is normalized such that the 
maximum value becomes 1. 

Step S3: Calculate the square of the variable 

Xi. 

20 Step S4: Set d=6 and substitute the constant 

K(d) pick up from t = 6 on the cosine side of Fig. 
6 for Q(6) . 

Step S5: Calculate d=d- 2 . 

Step S6: Check whether d is not less than 0. 
25 If It is not less than 0, proceed to step S7, and 
if it Is minus, proceed to step S8 . 

Step S7: Sequentially calculate the cyclic 
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equation for d=4, 2, 0, through the route process 
of steps S5 to S7, to sequentially find Q(4) , Q(2) 
and Q( 0 ) . 

Step S8 : This is the case that arithmetic of 
5 the cyclic equation is terminated because d 

becomes minus, and in this case, Q(0) which is 
found at the last in steps S5 to S7 is considered 
as cos 6 i . 

Step S9 : Output the value of cos 6 i found by 

10 the arithmetic process. 

Fig. 10 is a block diagram of a function 
arithmetic circuit of the present invention 
achieved by hardware, and the function arithmetic 
of the present invention achieved by the program 

15 control by the computer, which is shown by the 
embodiment of Fig. 5, is achieved by hardware . In 
Fig. 10, the function arithmetic apparatus of the 
present invention consists of the angle detection 
unit 10 and any apparatus 18, and the apparatus 

20 is provided with a function arithmetic circuit 19 . 
The function arithmetic circuit 19 consists of an 
arithmetic control unit 20, a constant table 22, 
a conversion adjustment unit 24 and a cyclic 
equation arithmetic unit 25. In the constant 

25 table 22, the shift umbers Sd and the constants 
Kd of sine and cosine, which are needed for 
calculation up to, for example, 45 degrees In the 
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constant table 16 of F±g. 6, Is stored for 
necessary values out of d, for example, d=l, 3, 
5, 7 of s±n and d = 0, 2, 4, 6 of cos. The conversion 
adjustment unit converts tlie address 1 input from 
5 the angle detection unit 10 Into the variable Xi, 
adjusts the shift number such that a maximum value 
of the variable Xi becomes 1 within a transform 
range and outputs it* The cyclic arithmetic unit 
25 calculates the cyclic equation obtained by 

10 transforming the Taylor series equation which 
calculates the sine function and the cosine 
function, in which the new known number Qd is 
considered as the number obtained by multiplying 
the known number d by the square of the variable 

15 Xi and adding the constant Kd after shifting the 
results by the shift number Sd. The arithmetic 
control unit 20 forces the cyclic equation 
arithmetic unit 25 to sequentially perform cyclic 
arithmetic based on the selection of the constant 

20 Kd and the shift number Sd of the constant table 
22, from higher order term to lower order term, 
for the number of terms of the Taylor series 
equation defined in advance when the variable Xi 
is output from the conversion adjustment unit 24 

25 and forces it to calculate the sine function sin 
0 i and the cosine function cos 9 i of the address 
i. 
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Fig. 11 ±s a circuit block diagram of a sine 
function arithmetic circuit used in Fig. 10. In 
Fig. 11, the sine function arithmetic circuit 25-1 
consists of a conversion shift circuit 26, 
5 multipliers 27, 30 and 40, an adder 34 and 

selection circuits 35, 36 and 38, and a shift 
circuit 32 and the adder 34 constitute a cyclic 
processing unit 28 which repeats calculation of 
the cyclic equation. The sine function 

10 arithmetic circuit 25-1 is provided with an 

arithmetic control unit 20-1. The conversion 
shift circuit 26 converts the input address i into 
the variable of Xi=i/(N/8) by shifting the bits 
by (N/8). The multiplier 27 finds the variable 

15 Xi 2 . The cyclic processing unit 28 performs the 
multiplication of Xi 2 and Q(d+2) in the cyclic 
equation of the equation (1). The shift circuit 
32 performs the shift of the shift number S(d+2) 
in the equation (1). The adder 34 performs the 

20 addition for adding the result of the shift 

circuit 32 to the constant Kd. The output of the 
adder 34 returns to the selection circuit 35. The 
selection circuit 35 selects the constant K7 in 
the first multiplication in the multiplier 30, 

25 selects the known number Q5 output from the adder 
34 in the next multiplication and selects the 
known number Q3 output from the adder 34 in the 
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third multiplication. The selection circuit 36 
selects the shift number S7 = 2" 6 in the first shift, 
selects S7=2" 5 in the second shift operation and 
selects S7=2" 3 in the third shift operation. The 
5 selection circuit 38 selects the constant K5 in 
the first addition, selects the constant K3 in the 
second addition and selects the constant Kl in the 
third addition. The multiplier 40 multiplies the 
known number Ql, which is output after three (3) 

10 cyclic processes for d=5, 3, 1 in the cyclic 

processing unit 28 are completed, by the variable 
Xi, which is output by the conversion shift 
circuit 26, and outputs the result of the 
multiplication as sin 0 1. This sine function 

15 arithmetic circuit 25-1 of the present invention 
only needs arithmetic processes of five (5) 
multiplications and three (3) additions by 
repeatedly using the cyclic processing unit 28. 
Contrary to this, in the case of the conventional 

20 example shown in Fig. 4, eight (8) multiplications 
and three (3) additions are needed, and the sine 
function arithmetic circuit 25-1 using the cyclic 
equation according to the present invention can 
find the result with fewer arithmetic processes. 

25 Also, since the selection circuit and the shift 
circuit has smaller circuit scale as hardware than 
the addition circuit and the multiplication 
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circuit, ±n the sine function arithmetic circuit 
25-1 according to the present invention, the 
hardware scale can be configured smaller and 
easier than the conventional hardware . 

Fig. 12 is a circuit block diagram of a cosine 
function arithmetic circuit 25-2 in Fig. 10. The 
cosine function arithmetic circuit 25-2 consists 
of a conversion shift circuit 45, multipliers 46 
and 50, an adder 46, a shift circuit 52, an adder 
5 4 and selection circuits 5 5 , 5 6 and 5 8 - Among 
these, the multiplier 50, the shift circuit 52 and 
an adder 54 constitute the cyclic processing unit 
48. Also, an arithmetic control unit 20-2 is 
provided, which controls an arithmetic process of 
the cosine function arithmetic circuit 25- 2 . 
When the address i is input, the cosine function 
arithmetic circuit 25-2 converts it into the 
variable Xi=i/(N/8) with a (N/8) shift operation 
by the conversion shift circuit 45, and then the 
variable XI 2 is found by the multiplier 46. The 
output of the multiplier 4 6 is input to the 
multiplier 50 of the cyclic processing unit 48 and 
multiplied by the known number Q(d+2) from the 
selection circuit 55. The output of the 
multiplier 50 is shifted in the shift circuit 52 
based on the shift number S ( d+2 ) from the 
selection circuit 56. The output of the shift 



circuit 52 ±s added in the adder 54 to the constant 
Kd selected from the selection circuit 58. The 
output of the adder 54 is input to the other 
selection circuit 55 and becomes final arithmetic 
5 output. The arithmetic control unit 20-2 

performs cyclic arithmetic such that d=6, 4, 2, 
0 is achieved following cyclic arithmetic of the 
equation (1) in the cyclic processing unit 48. In 
other words, in the first addition in the adder 

10 50, the selection circuit 55. selects the constant 
K6 corresponding to d=6 to perform the addition, 
and subsequently, the selection circuit 56 
selects the shift number S6=2-6 for the shift 
circuit 52 to perform the shift, and the selection 

15 circuit 58 selects the constant K4 for the adder 
54 to perform the addition output of the known 
number Q4 . This addition output Q4 is returned 
and input to the selection circuit 55, selected 
for the second multiplication in the multiplier 

20 50 and multiplied by Xi 2 . Subsequently, the 
selection circuit 56 sets S4 = 2" 4 to the shift 
circuit 52 to perform the shift process, and the 
known number Q2 is found in the adder 54 by the 
addition with the constant K2 selected in the 

25 selection circuit 58. This addition output Q2 is 
input to the selection circuit 55 again, and in 
the third addition, the multiplier 50 multiplies 
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XI 2 by the known number Q2 from the selection 
circuit 55, and after the shift of S2 = 2 _1 selected 
by the selection circuit 56 is performed in the 
shift circuit 52 # the addition with the constant 
5 KO which is currently selected by the selection 
circuit 58 is performed in the adder 54 to find 
the addition output QO . Since this addition 
output QO is a cosine function which is finally 
calculated, this is output as the arithmetic 

10 result. 

Fig. 13 is a flowchart of another embodiment 
of the sine function arithmetic process of the 
present invention in accordance with the function 
arithmetic processing unit 12 of Fig. 5 which 

15 keeps arithmetic accuracy by rounding arithmetic 
results to which protect bits are added. In this 
sine function arithmetic process of Fig. 13, for 
the known number Qd, the variable XI 2 , an 
Intermediate value after the shift and the 

20 constant Kd, the cyclic calculation is performed 
with the number of bits which is obtained by adding 
the number of protect bits (guard bits) to the 
final number of the bits, and then the necessary 
final number of bits can be ensured by rounding 

25 the protect bit portion of the final value. With 
these arithmetic processes in which the protect 
bit portion of the final value is rounded after 
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the cyclic calculation is performed with adding 
the number of the protect bits, it is possible to 
perform each calculation with the same bit 
accuracy as fixed-point arithmetic, and in the 
software process of Fig. 13, the processing load 
can be reduced because it is fixed-point 
arithmetic, and the circuit configuration in the 
case that it is achieved by hardware can be 
simplified. Assuming that bit data is b; protect 
bit data is g; the necessary final number of bits 
after the decimal point is eight (8) bits; and 
protect bits are two (2) bit, each of data in the 
sine function arithmetic process of Fig. 13 can 
be represented as follows. 

Constant Qd = sb.bbbb bbbb gg 

Variable XI = b.bbbb bbbb gg 

Square of variable XI = b.bbbb bbbb gg 
Constant Kd = sb.bbbb bbbb gg 

In the constant table 16 of Fig. 6, when the cosine 
function is calculated in the range from 0 to 90 
degrees, shift S2 is S2 = -l. The shift S2 = -l means 
shifting to left by one (1) bit, and in order to 
ensure the same accuracy with the case of 
calculating from 0 to 45 degrees, the calculation 
must be performed with the protect bits (guard 
bits) which are increased one (1) bit from the case 
of calculating from 0 to 45 degrees. The process 
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procedure of the sine function arithmetic process 
of Fig. 13 is described as follows. 
Step SI: Input the address 1. 

Step S2: Find the variable XI by the shift 
5 process of the address 1. At this point, two (2) 
bits of the protect bits are added to the variable 
XI. 

Step S3: Square the variable XI and round the 
result with half adjust. 

10 Step S4 : Set d= 7 and define the known number 

Q(d) as the constant K(7) which is read from the 
constant table 16 of Fig. 6. 

Step S5: Calculate d=d- 2 . d is varied to 5, 
3, 1 and 0 by the cyclic process. 

15 Step S6 : Check whether d Is not less than 1. 

If it is not less than 1, proceed to step S7, and 
If less than 1, proceed to step S8. 

Step S7: Multiply the variable Xi 2 by the 
known number Q(b+2) and round the results with 

20 half adjust to obtain the Intermediate value W* 
Perform the shift process of the intermediate 
value W with the shift number S(d+2). Round the 
value obtained by adding the constant K(d) to the 
intermediate value W with half adjust to obtain 

25 new known number Q(d). This calculation is 
repeated three (3) times for d=5 , 3, 1. 

Step S8: This is the case that the cyclic 
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calculation is terminated because d becomes less 
than 1 # and sin 0 1 is found by multiplying Q(l) 
by the variable Xi, and the value of the final 
number of bits is found by rounding the results 
5 with half adjust. 

Step S9: Output sin 0 1 as the arithmetic 
result. 

Fig. 14 is a flowchart of the cosine function 
arithmetic process according to the present 

10 invention which keeps arithmetic accuracy by 

rounding arithmetic results to which protect bits 
are added, and the process procedure is as 
follows • 

Step SI: Input the address 1. 

15 Step S2: Find the variable XI by the shifting 

address I by (N/8). At this point, two (2) bits 
are added to the variable XI as the protect bits. 

Step S3: Find the square of the variable XI 2 
and round it with half adjust* 

20 Step S4 : Set d=6 and define the known number 

Q(6) as the constant K(6) which is found from the 
constant table 16 of Fig. 6. 

Step S5 : Calculate d=d-2. In this case, d 
is varied to d=4 , 2, 0. 

25 Step S6: Check whether d is not less than 0. 

If it is not less than 0, proceed to step S7, and 
if it is minus, proceed to step S8. 
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Step S7: Multiply the variable XI 2 by the 
known number Q(b+2) and find the intermediate 
value W by rounding with half adjust . Then, after 
performing the shift process of the intermediate 
value W with the shift number S(d+2), add the 
constant K(d) and perform rounding with half 
adjust to find new known number Q(d). This 
arithmetic process is executed three (3) times for 
Q(4) , Q(2) and Q(0) through the loop process from 
S5 to S7. 

Step S8 : This is the case that the cyclic 
arithmetic is terminated because d becomes minus, 
and cos 0 i with the final number of bits is found 
by rounding Q( 0 ) found in step S7 with half adjust . 

Step S9: Output cos 6 i as the arithmetic 
result . 

For rounding in the processes of Fig. 13 and Fig. 
14, the half adjust Is taken as an example, but 
cutoff of lower bits may also be used. For the 
half adjust, by adding 1 to the position which is 
one (1) bit lower than the lowest bit position to 
be found, calculating and then performing cutoff, 
the half adjust can be achieved. Cutoff may be 
achieved by leaving the necessary number of bits 
after the decimal point and cutting off subsequent 
bits. 

Fig. 15 is another embodiment of the function 
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arithmetic circuit of the present invention based 
on the hardware configuration shown in Fig. 10 and 
an embodiment of a sine function arithmetic 
circuit which performs pipeline arithmetic. 
5 This pipeline sine function arithmetic circuit 
keeps arithmetic accuracy by rounding the 
arithmetic results to which protect bits are added, 
as is the case with the sine function arithmetic 
process of Fig. 13. 

10 In Fig. 16, a sine function pipeline 

arithmetic circuit 25-3 consists of a conversion 
shift circuit 64, a multiplier 65, a first 
arithmetic unit 66, a second arithmetic unit 74, 
a third arithmetic unit 80, a multiplier 86 and 

15 a rounding circuit 88. The first arithmetic unit 
66, the second arithmetic unit 74 and the third 
arithmetic unit 80 are in the pipeline connection. 
The first arithmetic unit 66 is comprised of a 
multiplier 68, a shift circuit 70 and an adder 72; 

20 the second arithmetic unit 74 is comprised of a 
multiplier 75, a shift circuit 76 and an adder 78; 
the third arithmetic unit 80 is comprised of a 
multiplier 82, a shift circuit 84 and an adder 85; 
and each of these three (3) arithmetic units 66, 

25 74 and 80 is provided with the same circuit 

configuration. The multipliers 65, 68, 75, 82 
and 86 and the adders 72, 78 and 85 are provided 
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with a rounding function (Round) and perform 
cutoff or half adjust for the arithmetic results. 
The cutoff is performed by cutting off lower bits 
after the valid bits. In the case of the half 
5 adjust, by adding 1 to the position which is one 
(1) bit lower than the lowest bit position of the 
valid bit , calculating and then performing cutoff, 
the half adjust can be achieved. In this case, 
for example, by adding appropriate protect bits 

10 (guard bits) to eight (8) bits after the decimal 
point which are valid bits and performing 
calculation, for example, by adding two (2) guard 
bits to each of the known number Qd, the variable 
XI and Xi 2 and the constant Kd and performing 

15 calculation, the calculation can be performed 
with practical arithmetic accuracy even in the 
case that rounding is the cutoff process. As the 
shift circuit 70, 76 and 84, circuits do not 
particularly have to be provided in actual 

20 hardware, and when the output signals of the 

multipliers 68, 75 and 82 are supplied to each of 
the adder 72, 78 and 85 as the input signals, only 
the function for changing the bit position may be 
provided. In this sine function pipeline 

25 arithmetic circuit 25-3, when the address i 

indicating angles is sequentially input in sync 
with the arithmetic cycles in conformity with a 
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pipeline control unit 20-3 , pipeline arithmetic 
is achieved , in which , for example , the first 
arithmetic result of the sine function is output 
in the ninth cycle, and after that, the arithmetic 
5 results of the sine function are sequent! ally- 
output for each cycle* 

Fig. 16 shows an embodiment of a cosine 
function arithmetic circuit according to the 
present invention which performs pipeline 

10 arithmetic. This cosine function pipeline 

arithmetic circuit 25-4 consists of a conversion 
shift circuit 95, a multiplier 96, a first 
arithmetic unit 98, a second arithmetic unit 105 
and a third arithmetic unit 112, the first 

15 arithmetic unit 98, the second arithmetic unit 105 
and the third arithmetic unit 112 are in the 
pipeline connection. The first arithmetic unit 
98 consists of a multiplier 100, a shift circuit 
102 and an adder 104; the second arithmetic unit 

20 105 consists of a multiplier 106, a shift circuit 
108 and an adder 110; the third arithmetic unit 
112 consists of a multiplier 114, a shift circuit 
115 and an adder 116; and three (3) arithmetic 
units have the same circuit configuration and are 

25 in the pipeline connection. The multipliers 96, 
100, 106 and 114 and the adders 104, 110 and 116 
are provided with a function for performing the 
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cutoff or the half adjust to the necessary number 
of the bits by rounding. In this cosine function 
pipeline arithmetic circuit 25-3, when the 
address i as angle information is sequentially 
5 input in each arithmetic cycle by pipeline 

arithmetic control of a pipeline control unit 2 0 - 4 , 
the first arithmetic result of the cosine function 
can be output in the eighth cycle, and after that, 
the arithmetic results will be sequentially 

10 output for each cycle. 

Figs. 17A and 17B are block diagrams of the 
case that the pipeline sine function and cosine 
function arithmetic circuits shown in Fig. 15 and 
Fig. 16 are applied to a fast Fourier transform 

15 circuit whose radix is 2. In Figs. 17A and 17B, 
the fast Fourier transform circuit with radix 2 
has a four-stage structure of a sorting circuit 
120, a butterfly arithmetic unit 122, a sorting 
circuit 130, a butterfly arithmetic unit 132, a 

20 sorting circuit 134, a butterfly arithmetic unit 
136, a sorting circuit 138 and a butterfly 
arithmetic unit 140. In this embodiment, fast 
Fourier transform is performed by inputting 16 
points of sampling data, for example. For the 

25 time-series input signal X0(k) to the sorting 
circuit 120, the butterfly arithmetic is 
performed among data of X0(k+8) in which the 
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number k ±s separated by 8 In the butterfly 
arithmetic unit 122. In the next second-stage 
butterfly arithmetic unit 132, the butterfly 
arithmetic is performed among data Xl(k+4) in 
5 which the number k of the input data XI (k) is 
separated by 4 . Subsequently, in the third-stage 
butterfly arithmetic unit 136, the butterfly 
arithmetic is performed among data X2(k+2) in 
which the input data X2(k) is separated by 2. 

10 Finally, in the fourth-stage butterfly 

arithmetic unit 140, the butterfly arithmetic is 
performed among data X3(k+1) in which the input 
data X3(k) is separated by 1, and Fourier 
transform output X4(k) and X4(k+1) are output. 

15 Fig. 18 shows input and output data in each 

combined portion of the sorting circuits 120 , 130 , 
134 and 138 and the butterfly arithmetic units 122, 
132, 136 and 140 of the first stage to the fourth 
stage in the fast Fourier transform circuit of 

20 Figs. 17A and 17B. In Fig. 18, out of 16 points 
of sampling data, X0 ( 0 ) , X0(2), and X0(14) 

whose numbers are even and X0 ( 1 ) , X0 ( 3 ) , • * • # and 
X0(15) whose numbers are odd are respectively 
input to two (2) Inputs, and the butterfly 

25 arithmetic is performed among data separated by 
K=8, and the arithmetic data shown on the right 
side of the butterfly arithmetic unit 122 are 
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obtained. These data are sequentially input to 
the sorting circuit 130 on the second row as two 
(2) units of data, and with the butterfly 
arithmetic unit 132, the arithmetic data on the 
5 right side thereof is obtained. Subsequently, 
the output data on the second row are input to the 
third-stage sorting circuit on the third row; the 
arithmetic result on the right side of the 
butterfly arithmetic unit 136 is obtained; and 

10 after this result is input to the sorting circuit 
138 as the fourth-stage input on the fourth row, 
the result of the fast Fourier transform output 
shown on the right side is obtained by arithmetic 
with the butterfly arithmetic unit 140. 

15 Referring again to Figs. 17A and 17B, as shown, 

for example, in the first-stage butterfly 
arithmetic unit 122, the multiplier 124 performs 
twist coefficient multiplication which 
multiplies the input data by a twist coefficient 

20 W° obtained via a buffer 145. For this twist 

coefficient Increases, the number of types are 
Increased by 2n such as 1,2,4 and 8 in the order 
from the first stage to the fourth stage. Each 
of the twist coefficients on each stage is used 

25 eight (8) times on the first stage, four times on 
the second stage, twice on the third stage and once 
on the fourth stage for multiplication. 
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Multiplication of these twist coefficients and 
each of data on each stage of Fig. 18 is shown as 
follows . 

(The first stage) XO ( 8 ) -W°, X0(9) -W°, X0(10) • 
W°. X0(11)-W°, X0(12)'W°, X0(13)'W°, X0(14)'W°, 
X0( 15) 'W° 

(The second stage) XI ( 4 ) -W°, Xl(5) -W°, Xl(6) • 
W°, Xl(7)-W°, 

X1(12)*W 4 , X1(13)*W 4 , Xl(14)-W 4 , X1(15)'W 4 

(The third stage) X2(2)-W°, X2(3)'W°, 

X2(6)-W 4 , X2(7)*W 4 , 

X2(10)*W 2 , X2(11)*W 2 , 

X2(14)*W 6 , X2(15)"W 6 

(The fourth stage) X3(1)*W°, 

X3(3) -W 4 , 

X3(5) -W 2 , 

X3(7) -W 6 , 

X3(9) 'W 1 , 

X3( 11) -W 5 , 

X3(13) *W 3 , 

X3(15) -W 7 

These twist coefficients used in the butterfly 
arithmetic units 122, 132, 136 and 140 are 
calculated in each of function arithmetic units 
144-1 and 144-2 in real time, based on the address 
1 given from address calculation units 142-1 and 
142-2. Each of the twist coefficient arithmetic 
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units 144-1 and 144-2 ±s provided with the sine 
function pipeline arithmetic circuit 25-3 shown 
in Fig. 15 and the cosine function pipeline 
arithmetic circuit 25-4 shown in Fig. 16, 
5 calculates the necessary twist coefficients in 
real time in sync with the arithmetic cycles of 
the butterfly arithmetic units 122, 132, 136 and 
140 on four (4) stages and stores necessary types 
of the twist coefficients into each of the buffers 

10 145, 146, 148 and 150, and on the fourth stage, 
eight (8) types of the twist coefficients W° , 
W 4 , and W 7 are sequentially calculated by the 

dedicated twist coefficient arithmetic unit 
144-2. When represented using W p , the twist 

15 coefficient is represented by multiplication of 
complex numbers as following equation. 

W P = e -J 2 *P/ H . oos(27Cp/N) - j • sin( 2 7C p/N) 
In this multiplication of complex numbers, the 
cosine function and the sine function must be 

20 calculated with an address p as an Input, and this 
is achieved by the pipeline cosine function 
arithmetic circuit and sine function arithmetic 
circuit of the present invention. The twist 
coefficient arithmetic unit 144-1 calculates the 

25 necessary twist coefficients by the arithmetic 
timing of the butterfly arithmetic of each stage 
in the fast Fourier transform and stores these 
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into the buffer 145, 146, 148 and 150. In other 
words, the twist coefficient arithmetic unit 
144-1 calculates W° , W° , W 4 , W° , W 2 , W 4 and W 6 as 
a twist coefficient sequence 300 as time t elapses 
5 and stores the twist coefficients indicated by 
black borders into the buffer registers 145, 146 
and 148 such that one (1) type is stored for the 
first stage, two (2) types for the second stage 
and four (4) types for the fourth stage. The 

10 twist coefficient arithmetic unit 144-2 is 
provided exclusively for the fourth-stage 
butterfly arithmetic unit 140, sequentially 
calculates W°, W 4 , W 2 , W 6 , W 1 , W 5 , W 3 and W 7 as a 
twist coefficient sequence 302 as time t elapses 

15 and stores this into the buffer 150. Then, the 
butterfly multiplication will be performed, in 
which the underside input data of sets of two (2) 
input data shown on the first stage to the fourth 
stage of Fig. 18 are multiplied by the 

20 corresponding twist coefficients. 

Figs. 19A and 19B are block diagrams in which 
the fast Fourier transform circuit of Figs. 17A 
and 17B is in a parallel configuration. In this 
fast Fourier transform circuit, the radix is 2 and 

25 the sampling points are 16 points, as is the case 
with Figs. 17A and 17B. In Figs. 19A and 19B, the 
fast Fourier transform circuit which is in a 
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parallel configuration consists of a first-stage 
butterfly arithmetic unit 152, a second-stage 
butterfly arithmetic unit 154 , a third-stage 
butterfly arithmetic unit 156 and a fourth-stage 
5 butterfly arithmetic unit 158. The first-stage 
butterfly arithmetic unit 152 consists of a twist 
coefficient multiplication unit 160 , comprising 
eight (8) multipliers, and a 16-plece adder group 
162. The second-stage butterfly arithmetic unit 

10 154 consists of a twist coefficient 

multiplication unit 164, which uses multipliers 
divided into two (2) sets of four (4) pieces, and 
an adder group 166, comprising 16 adders. The 
third-stage butterfly arithmetic unit 156 

15 consists of a twist coefficient multiplication 
unit 168, comprising four (4) sets of two (2) 
multipliers, and an adder group 170, comprising 
16 adders. The fourth-stage butterfly 
arithmetic unit 158 consists of a twist 

20 coefficient multiplication unit 172, comprising 
eight (8) multipliers, and an adder group 174, 
comprising 16 adders. In this fast Fourier 
transform circuit in a parallel configuration, 
the number of arithmetic stages from an input 

25 stage to an output stage is eight (8) stages, and 
the arithmetic process can be completed in eight 
(8) clock cycles. 
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A twist coefficient arithmetic circuit of Fig. 
20 is provided, corresponding to the fast Fourier 
transform circuit in a parallel configuration of 
Figs. 19A and 19B, in order to supply the necessary 
5 twist coefficients in real time to each of the 
twist coefficient multiplication units 160, 164, 
168 and 172 on each stage. The twist coefficient 
arithmetic circuit of Fig. 20 is basically the 
same as the fast Fourier transform circuit of Figs. 

10 17A and 17B. Corresponding to the first-stage 
butterfly arithmetic unit 152, the second-stage 
butterfly arithmetic unit 154 and the third-stage 
butterfly arithmetic unit 156, an address 
calculation unit 142-1 and a twist coefficient 

15 arithmetic unit 144-1 are provided. 

Corresponding to the first-stage butterfly 
arithmetic unit 152, the second-stage butterfly 
arithmetic unit 154 and the third-stage 
arithmetic unit 156 , buffer registers 176 , 178 and 

20 180 are also provided, respectively. For the 

fourth-stage butterfly arithmetic unit 158, an 
Independent address calculation unit 142-2 and a 
twist coefficient arithmetic unit 144-2 are 
provided, and a buffer 182 is also provided. The 

25 twist coefficient arithmetic unit 144-1 outputs 
the twist coefficients W°, W 4 , W 2 and W 6 at the 
timings of times to, tl, t2 and t3 indicated in 
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an output series 184 based on t = 0, 4, 2, 6 from 
the address calculation unit 142-1, stores W° into 
the buffer 176, stores W° and W 4 into the buffer 
178 and stores the twist coefficients W°, W 4 , W 2 
and W 6 into the buffer 180. On the other hand, 
the twist coefficient arithmetic unit 144-2 
sequentially calculates the twist coefficients W° , 
W 4 , W 2 , W 6 , W 1 , W 5 , W 3 and W 7 at the timings of times 
tO to t7 as shown in an output series 186 based 
on addresses t=0, 4, 2, 6, 1, 5, 3, 7 from the 
address calculation unit 142-2 and stores these 
into the buffer 182. In the fast Fourier 
transform circuit in a parallel configuration of 
Figs. 19A and 19B, eight multipliers provided in 
the twist coefficient multiplication unit 160 of 
the first-stage butterfly arithmetic unit 152 
performs the twist coefficient multiplication 
using the twist coefficient WO in the buffer 
register 176 at the timing of the time tl. The 
twist coefficient multiplication unit 164 of the 
second-stage butterfly arithmetic unit 154 
performs the twist coefficient multiplication 
using the twist coefficients WO and W4 stored in 
the buffer register 178 at the timing of the time 
t3. The twist coefficient multiplication unit 
168 of the third-stage butterfly arithmetic unit 
156 performs the twist coefficient 
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multiplication using the twist coefficients WO , 
W4 , W2 and W6 already stored in the buffer register 
180 at the timing of the time t5 . The 
fourth- st age butterfly arithmetic unit 158 
performs the twist coefficient multiplication 
using eight (8) types of the twist coefficients 

which are the twist coefficients W° , W 4 , and 

W 7 stored in the buffer register 182 at the timing 
of the time t7. In this way, by combining the fast 
Fourier transform circuit in a parallel 
configuration of Figs. 19A and 19B with the twist 
coefficient arithmetic circuit of Fig • 2 0 , the 
fast Fourier transform output can be obtained by 
performing the butterfly arithmetic across four 
( 4 ) stages in parallel in eight ( 8 ) arithmetic 
cycles from the Input for 16 points of sampling 
data of variables X0(0) to X0(15). In the 
arithmetic circuits of the sine function and the 
cosine function of Fig. 11 and Fig. 12, one (1) 
cyclic processing unit is provided and used 
repeatedly to the perform function arithmetic, 
but it is possible to convert this into the 
pipeline type as shown in Fig. 15 and Fig. 16. It 
is also possible to convert the pipeline sine and 
cosine arithmetic circuits shown in Fig. 15 and 
Fig. 16 into the function arithmetic circuits as 
shown in Fig. 11 and Fig. 12, in which one (1) 
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cyclic processing unit is provided and used 
repeatedly. 

In above embodiments, the sine function and the 
cosine function are taken as an example of the 
5 arithmetic process which is processed by 

transforming the Taylor series equation into the 
cyclic equation, but it is possible to apply to 
other trigonometric functions, such as a tangent 
function, a cosecant function, a secant function 

10 and cotangent function , and a logarithmic 

function and an exponential function included in 
general transcendental functions other than 
trigonometric functions to which the Taylor 
series equation is applied. When the Taylor 

15 series equation for the general transcendental 
functions other than trigonometric functions are 
transformed into the cyclic type of the present 
invention, the variable may be simply defined as 
X, rather than X 2 , and the Taylor series equation 

20 may be transformed into the cyclic equation which 
multiplies the known number Q and the variable X, 
shifts the result by the shift number S and then 
adds variable K to obtain new known number Q. 

The present invention has been described such 

25 that the hardware circuit can be easily achieved 
by defining the variable XI such that its maximum 
value becomes 1 and the constant Kd to be not 
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greater then 1, but considering the configuration 
of the arithmetic processing circuit or movement 
of the decimal point position, it is not limited 
to this value. In the case of hardware using the 
floating-point arithmetic circuit or the case of 
the software process by CPU or DSP which can use 
the floating-point arithmetic, by only 
transforming the Taylor series equation into the 
cyclic equation, it is possible to cut down the 
number of times of arithmetic, to reduce the 
processing load and to achieve the high-speed 
processing. 

Although in the above embodiments the case 
has been taken as an example where the arithmetic 
results of the sine function and the cosine 
function are applied to the fast Fourier transform, 
the present invention is not limited thereto but 
is directly applicable to any suitable 
apparatuses which utilize the arithmetic results 
of the sine function and the cosine function. 



